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A generalized Anderson single-impurity model with off-site Coulomb inter- 
' actions is derived from the extended three-band Hubbard model, originally 

proposed to describe the physics of the copper-oxides. Using the abelian 
bosonization technique and canonical transformations, an effective Hamilto- 
nian is derived in the strong coupling limit, which is essentially analogous to 



o 



^ ' the Toulouse limit of the ordinary Kondo problem. In this limit, the effective 

Hamiltonian can be exactly solved, with a mixed valence quantum critical 
point separating two different Fermi liquid phases, i.e. the Kondo phase and 
the empty orbital phase. In the mixed valence quantum critical regime, the 
local moment is only partially quenched and X-ray edge singularities are gen- 
erated. Around the quantum critical point, a new type of non-Fermi liquid 
behavior is predicted with an extra specific heat Gimp ~ T^/^ and a singular 
spin-susceptibility 'Ximp ^ 

same time, the effective Hamilto- 
nian under single occupancy is transformed into a resonant-level model, from 
which the correct Kondo physical properties (specific heat, spin susceptibility, 
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and an enhanced Wilson ratio) are easily rederived. Finally, a brief discus- 
sion is given to relate these theoretical results to observations in UPd^Cu^-x 
{x = 1, 1.5) alloys, which show single-impurity critical behavior consistent 

with our predictions. 
75.20.Hr, 72.10.Bg, 71.28.+d 
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I. INTRODUCTION 



Correlated electron systems ranging from heavy-fermions to high-temperature supercon- 
ductors exhibit a fascinating variety of anomalous behavior. An issue of considerable current 
debate is to what extent these systems in the normal state can be described as Fermi liquids 
(FL) where the low-energy excitations have a one-to-one correspondence to those of a nonin- 
teracting Fermi gas with the well-known behavior of a specific heat C = 7T, a temperature 
independent Pauli susceptibility and a low temperature electrical resistivity p ^ AT^. 
Most previous investigations on heavy-fermion systems, in particular, on the mixed valence 
problem are carried out within the framework of the FL theory. In the meantime a phe- 
nomenological marginal Fermi liquid (MFL) theory [||], proposed for high-temperature su- 
perconductors to explain certain features of the observed diversified anomalies, has brought 
about a lot of new developments in this active field of research. 

Recently, a lot of strong evidence of non-FL behavior has been reported for several 
heavy fermion systems 0, including UxYi^xPds and UPd^Cu^-x alloys [^-Q. In some 
of these systems, the non-FL behavior may be due to a T = quantum phase transition 
of cooperative origin p[, while in the others the singular behavior of thermodynamical 
and transport properties is consistent with a single-impurity critical scaling 0. Since no 
reliable methods have been devised so far to study the microscopic mechanism for the non- 
FL behavior in any lattice models, it is thus expected that the investigations of the related 
impurity models would shed some light on the basic physics of the lattice models, in the same 
sense as the original Anderson single-impurity model vs the single band Hubbard model (or 
Anderson lattice model). Various impurity models have been studied to explore the possible 
non-FL behavior [[7|-[T^. In particular, the two-channel Kondo system, intensively studied 
using different methods f7|-[T0| displays a kind of MFL behavior. However, it is not clear at 
all how to directly relate such a two-channel Kondo model to the real correlated f-electron 
heavy-fermion systems and high-temperature superconductors. In fact, a long standing 
puzzle of the mixed valence problem, namely to understand the complexity of phenomena 
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in realistic systems in terms of fundamental physical concepts [|T9| , is still not fully resolved. 
However, it seems to us that the key issue is whether the FL picture is correct for the mixed 
valence state. 



The mixed valence phenomena usually occur in rare-earth compounds in which the 
proximity of the 4f-level to the Fermi surface leads to substantial charge (valence) fluctu- 
ations and instabilities of the magnetic moment. Basically, the mixed valence state can 
be thought of as a mixture of two bonding states 4/"'(5(i6s)'" and 4/"~^(5(i6s)™"^^, which 
are nearly degenerate from the quantum mechanical point of view. This situation must 
be understood in terms of the hybridization of the above two configurations so that the 
hybridized level is only partially occupied in the low-energy excitations. Although such a 
situation bears a strong resemblance to the crossover of the single-impurity problem from 
free moment behavior at high temperature to a strong-coupling FL regime in the Kondo 
ground state. However, that crossover does not coincide with significant valence change in 
most mixed valence systems [^]. It was argued that the Friedel sum rule (in the sense of 
local charge neutrality) is important for the mixed valence problem and it should be satisfied 
for both valence states. In this sense the usual Anderson model is not complete (because 
the Friedel sum rule cannot be satisfied for two different valence states by varying only one 
parameter (the /— level)), so one should extend it to include the screening mechanism from 



the 5dQs electrons to describe the 4/ charge fluctuations [pi| , pj|] . 

Recently, the interest to this problem has been revived in connection with high Tc su- 
perconductors. Varma and others [p!^ -[T6| have considered a generalized Anderson model 
with finite-range Coulomb interactions, equivalent to a single-site version of the extended 
three-band Hubbard model proposed originally to describe the physics of copper-oxide su- 
perconductors. The numerical evidence for a mixed-valence quantum critical point came 
from a Wilson renormalization group study |T^, while a strong coupling limit was argued 



and considered for this critical point later [|T^]. Inspired by these studies, two of us ||n 
have shown that using the renormalization-group analysis a phase diagram with a mixed 
valence quantum critical regime separating the Kondo phase and the empty orbital phase. 



can be deduced. It has also been shown |jl^ that the strong couphng hmit can be derived 
straightforwardly from the original Hamiltonian using the canonical transformations and an 
exact solution to this model can be found in this limit. However, the physics behind the 
mixed valence state has not been fully revealed. The two basic important questions are: 
How to derive a complete effective model Hamiltonian in the strong coupling limit and what 
are the generic physical properties of the mixed valence state? 

In the present paper, we address these unsolved problems in the infinite-?7 generalized 
Anderson single-impurity model with finite-range Coulomb interactions. Through two suc- 
cessive canonical transformations in the bosonic representation, we can derive an effective 
Hamiltonian in a particular strong coupling limit, which is analogous to the Toulouse limit 
of the ordinary Kondo problem [|22|. In this limit, the effective Hamiltonian can be exactly 
solved, with a mixed valence quantum critical point separating two different Fermi liquid 
phases, i.e. the Kondo phase and the empty orbital phase. In this mixed valence quantum 
critical regime, the local moment is only partially quenched and X-ray edge singularities are 
generated. Around this critical point, a new type of non- Fermi liquid behavior is predicted 
with an extra specific heat Cimp ~ T^^^ and a singular spin-susceptibility Ximp ~ T~^/^. 
In the meantime, the effective Hamiltonian under single occupancy is transformed into a 



resonant-level model from which the correct Kondo physical properties (specific heat, 
spin susceptibility, and enhanced Wilson ratio) can be easily rederived. Finally, a brief 
discussion is given to relate our results with the experimental observations in UPd^Cu^^x. 
(x = 1, 1.5) alloys. 

The arrangement of the paper is as follows: We first derive the model Hamiltonian 
from the extended three-band Hubbard model in Section II. Then, in Section HI a non- 
perturbative approach for this single impurity model is formulated and a complete effective 
Hamiltonian is presented. In Section IV, the exact calculations of the physical properties 
in the strong-coupling limit are presented. Finally, discussions on experiments are given in 
Section V, while concluding remarks are made in Section VI. 
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II. DERIVATION OF THE LOCAL IMPURITY MODEL FROM THE 
EXTENDED THREE-BAND HUBBARD MODEL 



The extended three-band Hubbard model is based on the observation that the chemistry 
of the copper-oxides is unique in that charge fluctuations between the cations and the anions 
controlled by their mutual interactions occur at low-energies [^. Such dynamical effects 



were expected to qualitatively change the low-energy response, compared with the single 
band Hubbard model and its large Hubbard U limit, the t — J model. The model is composed 
of d hole states on the copper sites (one hole per copper in the undoped material) in (i^.2_^2 
orbitals which are strongly hybridized with p cr-orbitals on the neighboring oxygen sites. 
As the undoped material is a magnetic insulator, the holes on the copper sites must be 
localized by a reasonably strong on-site Coulomb interaction. Also included are the finite- 
range Coulomb repulsions between oxygen electrons on nearest-neighbor sites of the localized 
copper electrons. The corresponding Hamiltonian is given by: 

+ E (KkP\/i,<r + h.C.) + ^ E \k,k'ip\„Pk',a - Pk',aP\jiY. Aa'd^,-' " \) ^ (1) 

i,k,a i,k,k',a 

where pj: ^ is the Bloch representation of the oxygen orbitals, while V- ^ is the Fourier 
transform of the Coulomb interactions of the i— site copper with neighboring oxygens, and 
is the total number of copper-oxygen cells. The rest of notation has obvious meaning. 
As pointed out in the Introduction, the above lattice model has so far proven intractable 
to give us reliable information on the low-energy physics of the system. However, one 
can regard strong coupling in the lattice models as a low- dimensional critical phenomenon 
involving long-time fluctuations at each localized site, but no critical fluctuations in space. 
The long-time fluctuations at each site are independent because of the local conservation 
laws when coherent effects and inter-impurity interactions are weak. Therefore, one can 
use a single-impurity model to describe the copper-oxides or heavy-fermion metals under 
the physical hypothesis that the essential new physics is local in real space, on the scale 
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of a lattice constant, so all the interesting effects are in the time or frequency domain. If 
the singularities found in the impurity problem do not depend on any special symmetries 
lost when going from the single impurity to the lattice, they are likely to be relevant to 
the behavior of the full lattice problem as well. Consider now the single-site (copper site) 
version of Eq. (|ip in the following form: 

" fc,(7 

Y.{hv\/. + h.c.) + v^^pIp^^'. - ppA,J(T. dU' - 1)^ (2) 

where we set the single impurity at the origin of the coordinates and, as usual, assume 
the dx2-y2 orbitals are non-degenerate so that the quantum numbers {Id, frid) describing the 
orbital angular momentum, are fixed. It should be pointed out that this single-impurity 
model is still a three-dimensional system with two different angular momentum orbitals of 
the d-electrons and p-electrons. In the following, we have to reduce this three dimensional 
model to a one-dimensional model by projecting the p cr-orbitals onto the dx^-y^ orbitals. 

First of all, we transform the plane wave representation of the p-electrons to the spherical 
representation in view of the assumed spherical symmetry about the impurity site. 

Pk,a = J2<^\^^^> Ck,a,l,m, (3) 
l,m 

— * A 

where we have used the definition k = (k, k) and the quantum number m can take values 
— /, — / -|- 1, +1. If we project the p a-orbitals onto the d orbitals, the angular quantum 
number I is required to take only one value, that is the angular quantum number of the d 
orbital / = Id- In what follows, this quantum number is dropped, and the quantum number m 
is regarded as a " channel" index of the p orbital electrons. Since the dispersion relation of the 
p-electrons also has spherical symmetry, the kinetic energy term in Eq. cTi^k~f^)p\ Pk a 

is transformed into Z]fc,o-,m(^fc ~ f^)^^! am^k,a,m, which implies there are several conduction 
electron channels in the reduced model. At the same time, the hybridization term between 
the d- and p-electrons ^ llk,Shp\„da + ^-C-) is changed into Y.k,a{'^kCl^a,mado + h.c). 
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As required by symmetry, the localized d orbitals can hybridize only with p orbitals with 
the same angular quantum numbers. Due to the non-degeneracy of the d orbitals, there is 
only one channel with m = rrid of p-electrons hybridizing with the local d orbitals, which in 
the following is called hybridizing channel and the others are called screening channels. 

Second, the finite range Coulomb interaction between the localized orbitals and the 
delocalized orbitals, or the X-ray edge (XRE) like potential scattering, is considered. Due 
to the spherical symmetry about the impurity site, the scattering potential V^g depends 
only on the angle between k and k', so it can be expanded as 

'^Vk,k',i,m < k \ l,m >< l,m \ k' >, 

l,m 

where the only remaining term in the summation over the angular momentum is / = Id- 
Finally, putting everything together, the complete form of the reduced one-dimensional 
single-impurity Hamiltonian is obtained as: 

H = E (^fc - ^^)Cl^a,mCk,a,ra + E ^ddld„ + Ud\d'^d\dl + + ^.C.) 

+ 



1 1 

777 E E ^m{Cl^a^j3k',a,m- Ck',a,mCl,,^^)(^dl,da' - (4) 



k,k' ,a {m=—ld,...+lil} u' 

where, without loss of any generality, the hybridization strength t^, and the Coulomb in- 
teraction parameters Vk.k'.m have been further assumed to be momentum independent, as 
in the usual treatments of the Anderson single-impurity model. This is the so-called lo- 
cal copper-oxide model or generalized Anderson single-impurity model with finite-range 
Coulomb interactions used to describe the local version of non-FL properties of the copper- 



oxide compounds |T3 Hl5|] . Here we have presented an explicit derivation. Although the 
derivation itself is given in terms of copper-oxygen orbitals, valid only for cuprates, the 
charge fluctuation physics, incorporated here via including the rf— orbitals and finite range 
Coulomb interactions between d and p electrons, is also materialized in other mixed valence 
compounds, especially in heavy fermion systems. Therefore, this model can be applied to 
these systems as well. It is worth mentioning that the infinite dimensional technique has also 



been used to study the above extended three-band Hubbard model ||T^, and a single- channel 



generalized Anderson single-impurity model with self-consistent condition has been derived. 
Since the infinite-dimensional approximation overemphasizes the dimensional symmetries, 
it, probably, will miss some other important symmetry properties in the considered corre- 
lated electron systems, and the most important mixed valence physics we are concerned with 
(the screening effects) is not included in that generalized Anderson single-impurity model. 



III. BOSONIZATION FOR THE GENERALIZED ANDERSON 
SINGLE-IMPURITY MODEL 

Usually, the interesting low-energy physics contained in the above single- impurity model 
involves only two configurations: rid = and = 1, and we can remove the = 2 
configuration by letting [/ — oo. According to the numerical RG analysis for the usual 
Anderson impurity model, apart from the particle-hole symmetric case, reducing to the 
Kondo model, particular attention should be paid to the asymmetric case In the 

infinite-U limit, the model Hamiltonian is defined hj H = + Hg, where 

h h 

k,a a a ^ k a 



o- cr,(7' tJ-^f 

t 1 1 /l -j- 

fc,(7,m>0 cr,m>0 Z Z ^ ,^ ,^>g 

where C^-^m = -^J2kCk,a,m is the conduction electron operator at the origin of the coor- 
dinates, Ud = J2aT^d,a, N is the number of the lattice sites, and we have separated the 
Hamiltonian into the hybridizing (m = 0) and screening parts {m^O). The chemical po- 
tential of the conduction electrons has been chosen to be zero and the energy level of the 
local impurity is also assumed to be very close to it because of the character of the mixed 
valence state. We include a uniform external magnetic field h to calculate the total spin 
susceptibility enhancement. The infinite U limit keeps the low-energy physics of the model, 
but we should work under a local constraint < 1, which is really crucial for determining 
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the physical behavior of the low-energy excitations. Equally, we can rewrite Hh in another 
form, which is useful for the following analysis. 

+\/oE(Ci,oa,o - \){nd,. -\) + K,'E(^ioa,o - \){nd,, ~\) + V^Y.CUC,Ad.. (6) 

(T ^ ^ a ^ ^ a 

where Vq = + ^) is the parallel spin scattering potential, V^' = (V^ — x) is the opposite 
spin scattering potential, and V^ = ^ the spin-flip scattering potential. 

We first bosonize the screening part. Note that it has only one Fermi point for each 
channel and the dispersion is linearized ek = {k — kp)/p with a cutoff ko and p = {hvp)~^, 
so we can define the bosonic operators as p6|] . 



^ ku-k ku 
bk,a,m = rnrF E ^q,(T,m^q+k,a,m, ^\,a,m ~ /nTF E ^q,<J,m^q-k,c7,m, k > 0, 

VkN vfcA^g=fc 
which obey the standard commutation relations. Since the spin degrees of freedom of the 
screening channel electrons are trivially involved, they can be separated from Hg, which then 
becomes 

k l2k 1 

Hg = E -(4,m«fc,m + el^ek,m) + E ^"^V T7(4,m + ^fc,™) ("-d - -) 
k,m>0 P k,m>0 * ^ 

+ 1^ J^J dx[{ek,me^'^ + el^e-^'^n], (7) 

k>0,m ' 

where ak,m = -^{bk,i,m+bk,i,m) and 6^,^ = -^{bk,i,m-bk,i,m) are the charge- and spin-density 
operators, respectively. Moreover, we can assume Vm = K for all m > without loss of 
generality, so the channel index can be dropped. Thus, the screening part is simplified as a 
single spinless channel: 

H's = E -4«. + E Vs\l^ial + ak){nd - \ (8) 

where we have dropped the spin-density operators because they are not involved in the inter- 
actions with the local impurity, while the charge-density part is kept: = -^J2m>o(^k,m, 



Vg = \/2NsVs and Ns is the number of the screening channels. Using the inverse bosoniza- 
tion, we can transform bosons back to fermions: 
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Hi = ^ekslsk + ^YlisW - Sfe'4)K - (9) 



Interestingly enough, the strong electron-phonon interaction accompanying a valence change 
on a localized impurity bears the same form as H^. So in the abelian representation, Hs can 
be thought of including the effects due to the large cell volume change (15 percent) accom- 
panying the charge fluctuations, which is also expected to significantly alter the dynamics 
in the mixed valence problem. Therefore, the physics we will discuss in the strong-coupling 
limit to some extent reflects some essential aspects of the real mixed valence systems. 

The same bosonization procedures can be used for the hybridizing electrons as well, 
and the boson operators are defined as before. Restricted to the low-lying excitations, the 
hybridizing Hamiltonian is expressed as: 



k i h 

Hl= J2 -K,ah,a + earid + -{nd,^ - nd,i) + t ^{Cld^ + h.c.) + 14 E ClC^dld^ 

k>0,a P ^ 



fe>0,(T ' fe>0,<T ' 



Here Ca{x) = y/kr)exp{J2k>o ~^ibk,a^^''^ — ^kx^y jg fermion expression in terms of 



/kN 

the bosons. 

In order to derive an effective Hamiltonian in the Toulouse limit, we perform the following 
two canonical transformations 

C/ = exp{E^=K-4)K- S^exp{ E -^{h,a-bl ){nd,a-h}- 
k>oVkN 2 {k>o,a}ykN 2 

As a matter of fact, the first canonical transformation transfers the singularities of the 
screening channel into the hybridizing channel via the local impurity. The second canonical 
transformation not only makes the phase shifts due to hybridization and part of the parallel- 
spin XRE scattering compensate each other, but also simplifies the spin-flip scattering in 
the hybridizing channel. This way the hybridizing electrons become free for a special strong 
coupling limit, which shares the same physical meaning as the ordinary Toulouse limit of 
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the usual single-impurity single-channel Kondo model. The transformed model Hamiltonian 
is given by: 

~ k k 

H = -bl,ah,a + {ed-VQkD)nd + tY,{daSo + h.c.) + V±kDY,dld^ + Y.~^i"'k 

k>0,a P a a k>0 P 

^ fc>0,(T ' k>0,a ' 

+(K-^)E\/5(4 + «.)K-^), (11) 

P k>o * ^ 
Note that the |(?T-d,t — i^d,i) term in if/^ has been exactly canceled by the canonical trans- 
formation S, which is similar to the Emery-Kivelson approach for the two-channel Kondo 



system ||T^. A fermionic form is obtained by the inverse bosonization H^jj = Ht + SH, 
where 

h 

Ht = E(efe + ^(^)ClaCk,a + (ed - V^kD)nd + V^ko E 4da 

k,a o" 

+t'Y{d„so + h.c.) + E (^ks\sk, 

c k 

+(v;-^)(4so-^)K-^). (12) 

This effective Hamiltonian, when Vq = K = ^ and Vq = 0, reduces to Ht, which is the same 
as the strong coupling Hamiltonian argued by the renormalization-group analysis [|T^. It 
is obvious that Ht is analogous to the Toulouse limit of the ordinary Kondo problem. It 
should be emphasized that the full effective model Hamiltonian in the strong coupling limit 
is derived here exactly, from which the corresponding physical properties in different phases 
can be calculated explicitly. 

IV. PHYSICAL PROPERTIES OF THE MODEL HAMILTONIAN IN THE 
TOULOUSE STRONG COUPLING LIMIT 
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A. Exact solution of the Toulouse limit Hamiltonian 



Now, let us first study the strong-coupling limit Hamiltonian Ht with Vq = = ^ and 
Vq = 0, where only part of degrees of freedom of the local impurity a = :^('^t + ^l) 
coupled to the conduction electrons, while the remaining part P = ^{d^ — d^) is decoupled 



except for the constraint |14|, |15]. Thus, the effective Hamiltonian is: 



k,cr k 

+ {ed + V^kD)n^ + V2t{sla^ + h.c.), (13) 

where the local constraint Ua + np < 1 will play a crucial role in determining the behavior of 
the system, and thus we must handle it exactly. In impurity zero-occupancy, the hybridizing 
and screening electrons are completely free with a zero phase shift due to the Friedel sum rule. 
In impurity single-occupancy, the localized impurity disappears, but the total energy of the 
system has a shift (e^ — Vj_A;d), and according to the Friedel sum rule the hybridizing electrons 
and the screening electrons have a unitarity phase shift ^ per spin. This corresponds to the 
behavior of the Kondo model exactly at its PL fixed point. If we want to calculate the 
physical properties, the leading irrelevant interactions should be involved, equivalent to 
driving the system away from its fixed point 



When the impurity occupancy fiuctuates between zero and one, it requires that ~ 
—Vxkoi corresponding to the mixed valence quantum critical point found in Ref. [^]. The 
crucial point is that Ht conserves |T^|T^, which makes the following calculations avail- 



able. The strong coupling limit Hamiltonian in the subspaces of n^j = 0, 1 are the following, 
without any constraints 

Hn^=o = + + E ^kslsk + (ed + V^koW + ^t{sW + h.c), (14) 

k,(T k 

Hn,=i = E(efc + ^^)CUCk,a + E ^ksisk + {ed - V^ko). (15) 

k,a k 

Thus, we can carry out our calculations in the subspaces of n/j. Any eigenstates of Ht 
should be decomposed of the known eigenstates in the subspaces = 0, 1: 
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I ^ > = | >„^=o + I >np=l ■ (16) 

In the np = subspace, the constraint is satisfied for any a, so the quadratic Hamiltonian 
can be exactly solved, and a useful expectation associated with a is derived as: 

where = (ed + Kl/cd), F = pt^ is the resonance width of the a particle with the conduction 
electrons, and the phase shift of s-electrons at the Fermi point due to the presence of the 
local impurity a is found to be | — tan"^^^)- Generally, it is known that the impurity 
energy level is not exactly pinned at the chemical potential, being very close to it. However, 
for simplicity, we assume the extreme condition letting ec = 0, so the above phase shift is 
chosen 7r/2 approximately. In the subspace = 1, the condition = is required by the 
local constraint, so the conduction electrons become free. More importantly, the essential 
physics of this mixed valence quantum critical regime is contained in the Green function of 
(3 in the restricted Hilbert space, which can be expressed as 

= -e{r) < I /3(r)/3t(0) | >„^=o +0{-r) < | /3H0)/5(r) | >„^=i . (18) 

Due to the conservation of np, the first term can be regarded as an XRE absorption process, 
while the second term as an XRE emission process. From the well-known work on XRE 



problem [27], the retarded Green's function /5 is proved to be 



This result shows that the single-particle Green function of /3 impurity displays XRE singu- 
larities or Anderson orthogonality catastrophe at this mixed valence quantum critical point. 
Based on these results, we find that 

<r^^>~Tl/^ >=i, T -> 0. (20) 
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Thus, we have < >= 1/2, and from the Friedel sum rule, the local impurity approximately 
has a phase shift tt/A per spin at the Fermi point at zero temperature. Moreover, the 
correlation function of the charge density Pd = ""-a + — | can be evaluated as 

< Mr)MO) >« (^H-^f. r>0. (21) 

However, the correlation function of the spin density cr^ = q;^/5 + /?^q; is more subtle because 
of the local constraint. For r > 0, it is written as 



< ^dir)a'M >=< at(r)/3(r)/?t(0)a(0) >„^^„ + < /?t(r)a(r)at(o)/?(0) 



=1 • 



When the local constraint is carefully considered, the nonzero contribution to the first ex- 
pectation value is 

while the nonzero contribution to the second term is 

where Hq is the effective Hamiltonian in the = subspace, and Hi is the effective 
Hamiltonian in the rijs = 1 subspace. These expectation values can be evaluated in the same 
way as the retarded Green function of (3 p^. We thus derive the final result 



1 -jj-i^ 

< ^l{r)am >^ 2[- W\ T > 0. (22) 

1 sin^-ni T) 

The longitudinal spin-spin correlation function is dramatically different from the usual single- 
impurity model (1/t^ in the low temperature limit). Here we have assumed ~ 0. Since 
the local impurity decouples from the magnetic field and hybridizing electrons at this limit, 
the singular part of Ximp and Cjmp vanishes exactly. A perturbative calculation around the 
solvable strong coupling limit is needed to explore the generic behavior of the mixed valence 
phenomena. 
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B. Toulouse limit in the single-occupancy 



We now turn to the discussions away from the strong couphng hmit. As a test, we have 
to recover the well-known results of the single- impurity model in the Kondo limit. Since 
^d,T + = 1 in the Kondo Imit, we have the following effective Hamiltonian: 



k>0,a P c P k>0,a ' ^ 

fc>0,(T ' 

+^ E / ^^[(^M - KiY'^ + (^,T - KiY-''^ (23) 

where the hybridization term disappears due to the single occupancy constraint and we 
have dropped the screening channel because it decouples from the the impurity. Moreover, 
the charge degrees of freedom for the hybridizing electrons can also be separated out if we 
introduce the charge and spin-density operators as = :^(&fc,t+^fc,j,) and = ;^(^fc,T~^fc,j,)) 
respectively. After that, the effective Hamiltonian for the Kondo phase is obtained as: 

H^f = E -bih + V^ko E 4da + ^{Vo- K - ^) E + bk){na, - n.J 

k>oP <y V ^ P k>o^ 

+1^ E /I / dx[{h,me''^ + bl^e-^'^)]. (24) 
It is very useful to perform the following canonical transformation 

P = exp{-l Y: -^{h - bl){d\d^ - d\d^)}. (25) 
2 VkN 

When we transform the bosons back to fermions, -f^i// becomes 

+ - ^o' + El^lC-,, - C,,Ct)(n,,T - + ^(n,,^ - n.J. (26) 

Note that due to the single-occupancy constraint nd^^+nd^^ = 1, the local impurity is reduced 
to SU{2) spin-1/2 operators: 

d\d^ = S+, dld^ = S-, l(d\d^ - dld^) = S\ 
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In fact, the local impurity orbital can also be expressed as a U{1) spin-1/2 representations 
in terms of spinless fermions without any constraints: 

= /^ s- = /, s^ = ff - ^, 

which is in fact the one-site Jordan- Wigner transformation for a spin-1/2 Pauli operator. 
Thus, the final version of the Kondo effective Hamiltonian is 

+ [V2{V, - --) + ^]^E(Cfc - Ck'CDiff - \). (27) 
p p k,k' ^ 

This is a single-channel resonant-level model P3 , from which the well-known Toulouse limit 
Hamiltonian can be deduced when Vq — V^' = ^(1 — ^): 

H5f = E(^^ + ^)CiC, + V^^Y^iClf + fCu) + ^(/V - \)- (28) 

Here it is exactly derived from the Anderson single-impurity model with the single- 
occupation constraint. Using the Toulouse limit conditions, the Fourier transformation 
of the /-electron propagator Gfiuj) is evaluated as 

Gf{uj) = [uj-^ + iVKSign{uj)]~\ (29) 

where Vk = The corresponding density of states is 



1 



K 



This expression brings out the essential feature of the Kondo problem. In the strong coupling 
regime, the spin degrees of freedom of the local impurity have been quenched out by a 
conduction electron (forming a spin singlet), and other conduction electrons just hybridize 
with the impurity charge degrees of freedom with Lorentzian width Tx-, or they only feel a 
resonant scattering potential provided by the spin singlet at the origin. Thus, the change of 
the free energy due to this hybridization may be calculated by using the density of states: 

1 ^ K 

5Fimp = - dujf{uj)tarr^{ r-), (30) 
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where f{uj) = (e^"^ + 1) ^. Then, it is straightforward to evaluate the impurity contribution 
to the entropy limT^Qlimh^QSimp = 0, the specific heat Cimp = of— ^, and the spin suscep- 



tibihty Ximp 



27rr; 



-. Thus, the known universal Wilson ratio R 



w 



' TXtmy ^ l( TXbulk ' 



2 is 



recovered 28,29 



C. Generic properties in the mixed valence quantum critical regime 



Let us now concentrate on the physical properties of the mixed valence regime around 

■MV _ 
eff - 



the mixed valence quantum critical point. The effective Hamiltonian is given by Hf"^^ 



Hq + 5H, where 

Ho = Yui^k + -^'^)Ck,aCk,^ + ^ksisk + (cd - V±kD)ni3 

k,a k 

+ {ed + V^kD)na + V2t{sla^ + h.c.) 

Z Pa Z Z Z Pa 

+(v;--)(4so-^K (31) 

p ^ 

where ~ — VxA;^ and 5H is the perturbation away from the mixed valence quantum critical 
point. Obviously, there is a decoupling of charge and spin densities for the hybridizing 
electrons which makes our following calculations transparent. After that separation, we get 

h 

Ho = E(^fc + —7K)3k9k + E ^kPlvk + E efc4-5fc + - V^ko + V±kD)n^ 

k V2 

+ {ed - V^kD - V±kn)nfs + V2t(4«t + h.c.) 

5H = -^{Vo + - -KpIpo - l){p<i -\) + ^{Vo - K - -pMgo - l)^'d 

+ {Vs--){slso-hpd. (32) 
P ^ 

Here the spinless p— and ^f— electrons describe the charge and spin degrees of freedom of the 
hybridizing electrons, respectively. Hq corresponds to our mixed valence quantum critical 
point in the strong coupling limit. As follows from our calculated charge density and spin 
density correlation functions of the local impurity of Hq, only interactions associated with 



the spin density correlation function in the perturbations are singular. We also know the 
local propagator of the spin part of the hybridizing electrons 



G,{r) = - < Trgl{T)9o{Q) >= knTY. / dk, ^ ^ p- 



but ^3(0) ^ ^ for /i ^ 0. 

To the second order in perturbation, the singular contribution of SH to the impurity free 
energy reads as: 

1 /-VT 



= -77 / dT< 5H(r)5H(0) > 
z Jo 



= - Y{G^(0) l'^^ dr < al{T)am > + fj^ dTG%T) < al{T)am >}, (33) 

where '^=;^(^~^'~^)- It turns out that the first term yields the dominant contribution to 
Ximp-i while the second term yields the dominant contribution to Cimp- At low temperatures 
(T ^ r), the asymptotic forms of the propagators can be used in the expression of AFj^p 
to yield 

(34) 

All integrals must be regularized by introducing an ultraviolet cut-off F"^. Finally, the 
low-temperature singular behavior of the spin susceptibility and specific heat is found to be 

Atmp — p TtT ' '^"^P ~ 16 F ' 

where A = ^^^^^^^^^^^-^ and T(x) is the Euler Gamma function. Higher order perturbation 
terms in A yield subdominant contributions at low temperatures, and will not modify these 
results. Therefore, the behavior of susceptibility and specific heat obtained here is generic 
in the mixed valence regime. The spin susceptibility of the local impurity explicitly reveals 
the mixed valence physics around the strong coupling limit: not fully quenched magnetic 
moment and exhibition of the XRE singularities. Put it another way, the local moment 
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fluctuations are essential for the mixed valence state. When we consider the leading power- 
law singularity of the spin susceptibility and specific heat, we find the effective Wilson ratio 
Tximp/Cimp = 32/7r^, which should be observable in experiment. The present method allows 
us to see explicitly how the dominant singular behavior of the thermodynamic quantities 
is governed by the leading irrelevant interactions in the strong coupling Hamiltonian, and 
thereby correctly captures the generic properties of the generalized Anderson single-impurity 
model in the mixed valence phase. 

V. EXPERIMENTS ON UPd^Cu^^-j, ALLOYS 

Now we try to relate the above obtained results to real strongly correlated f-electron heavy 
fermion systems, e.g. U Pd^Cu^^x alloys in which U ions have partially-filled f-electron shells 
and give rise to magnetic moments interacting with the spin and charge of the conduction 
electrons. What is striking about these materials is the non-FL behavior of their physical 
properties which exhibit weak power law or logarithmic divergences in temperature and 
suggest the existence of a critical regime at T = 

The parent compound UCu^ is a prototype Kondo impurity system at high tempera- 
tures. As temperature is lowered, the intersite interactions drive a transition to a long-range 
antiferromagnetic (AFM) order of incompletely compensated uranum moments at IbK. But 
as Pd is substituted on the Cu sites, the AFM order is suppressed, with the ordering tem- 
perature being driven to zero. In particular, the recent neutron scattering experiments 
showed that the magnetic response function for x=l and x=1.5 has no appreciable momen- 
tum (except for the ionic form factor) and temperature dependence and has a scale invariant 
form S{lo) ~ uo~^/^ at very low energy scale (less than 10 mev). Moreover, the specific heat 
and spin susceptibility measurements also showed that between IK and lOK the temper- 
ature variations for x=l and x=1.5 Cimp/T ~ Ximp apparently deviate from the logarithmic 
temperature dependence (expected for a two-channel Kondo model at low temperatures), 
and can be much better represented by a power law Cimp/T ~ Ximp ~ T^"^ with A = 0.32. 
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Furthermore, the neutron data at various energies and temperatures have been fitted 
by a single scahng function, and it has been argued that the existence of such a scahng 
behavior provides a strong evidence in favor of single- impurity critical phenomena 0. The 
conformal dimension of the spin correlation extracted from this fitting A = 1/3 is consistent 
with the specific heat, magnetic susceptibility and resistivity data [Q. On the other hand, 
it is very clear that neither the single channel Kondo (predicting the FL behavior) nor the 
two-channel Kondo (leading to logarithmic singularities) can provide any acceptable expla- 
nation for such a striking behavior. It is thus natural to expect that the additional charge 
fluctuations inherent in the mixed valence phenomena, may shed some light on this puzzle. 

Here we propose a possible explanation for these experimental observations in terms of 
the generalized Anderson model considered in this paper. First, U Pd^Cu^^x with x=l and 
x=1.5 may be described as mixed valence states because the f-electrons of U-ion is only 
partially filled and the impurity energy level is split in the crystalline field into a large 
number of sublevels with possibly one of them being close to the chemical potential of the 
conduction electrons. Second, there are probably strong screening interactions from the 
off-site conduction electrons and strong impurity-phonon interaction due to the large lattice 
distortion accompanying the impurity occupancy changes. Finally, due to the presence 
of strong screening and phonon interactions, the motion of the U-ion must drag a heavy 
screening cloud with it. Hence its motion is very "slow", so the correlations from the 
inter-impurity interaction are rather weak. Thus, we can use our generalized Anderson 
single-impurity model to describe these materials. 

According to our calculations, we have found that the extra specific heat and spin suscep- 
tibility in the mixed valence quantum critical regime can be expressed as Cimp/T ~ Xmp ~ 
2"-3/4 ^YiQ strong-coupling limit. As for the magnetic response function ^(u;) from the 
neutral scattering measurements, we can obtain it after performing Fourier transformation 
for the spin-spin density correlation derived in the strong coupling limit of the mixed valence 
regime. We found that in the extremely low-energy scale S{uj) ~ uj^^^^. Apparently, the 
singularity exponent derived from our theoretical calculations (-3/4) is stronger than the 
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experimental results ~ —1/3. A possible way to recoincide these results is to note that in 
our treatment we have pushed the impurity scattering to the unitarity limit (phase shift 
7r/2), while in realistic systems the unitarity limit may not be reached, so the exponents 



might be different In view of the experimental relevance, it is very important to find 
other ways to treat this intriguing problem for cross-checking. 



VI. CONCLUDING REMARKS 

To summarize, in this paper we have considered a generalized Anderson model with 
screening channels, which can be considered as a single-site version of the three-band Hub- 
bard model originally proposed to describe the cuprates. Using the bosonization technique 
and canonical transformations we have found the strong coupling Toulouse limit of the ef- 
fective Hamiltonian which can be solved exactly. The physical properties around the mixed 
valence quantum critical point separating the two FL phases (the Kondo and the empty 
orbital phases) have been calculated explicitly to predict a new universality class of non-FL 
behavior. 

Our strong coupling treatment of the mixed valence problem provides the following 
physical picture for the mixed valence state. In the hybridizing channel, the phase shift due 
to hybridization is compensated by the parallel-spin XRE scattering via transformation S, 
while the XRE scattering in the screening channel is converted into an effective hybridization 
with the local impurity via transformation U. However, only "half of the local impurity 
{a particle) hybridizes with the screening electrons (in fact "holes") and gives rise to a 
standard FL behavior. The other "half of the local impurity {(3 particle) does not mix 
with the screening electrons and its expectation value < n/j >~ T^^^ at low temperatures. 
Nevertheless, its dynamic fluctuations lead to XRE singularities and a divergent power-law 
susceptibility which are consistent with observations in U PdxCu^-x alloys with x=l and 
x=1.5. 

The non-trivial physics shown in this paper is grasped due to a careful treatment of the 
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single occupancy constraint in the Toulouse limit. The essence of this constraint is to satisfy 
the unitarity condition, i.e., to use a complete set in the Hilbert space. Any mean- field type 
treatment of the constraint will certainly miss this basic point. 

We would hke to thank M. Aronson, M. Fabrizio, T. Giamarchi, A. C. Hewson, G. 
Kotliar, Y. L. Liu, and A.M. Tsvehk for helpful discussions. G.-M. Zhang is supported by 
the Science and Engineering Research Coucil (SERC) in Britain. 
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